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Introduction 


Consider an evolutionary PDE with one spatial variable 

Ut = F{u,Ux,. ■. 

and a symmetry of this equation, i.e. another system 

Us = G(u,Ux, . ■ ■jU^''^) 


( 0 . 1 ) 


( 0 . 2 ) 


commuting with the first one, 


{Us)t = {Ut)s- 


The set of the stationary points = 0 of the symmetry is a finite-dimensional invariant manifold for 
the system (0.1). Particularly, in important examples, the invariant manifold can be described as a set of 
stationary points of a first integral of the system (0.1): 


Su(x) 


II 

e 

dt 


In this case it is known ([BN], [Mo]) that the restriction of the initial PDE to the invariant submanifold is 
a Hamiltonian system of ODEs. In particular Bogoyavlenskii and Novikov [BN] found a universal scheme to 
construct the Hamiltonian function of the reduced system in terms of the Hamiltonian of the original PDE. 
In this paper we extend this scheme to more general finite dimensional invariant submanifolds specified 
by local x- and time-dependent symmetries and conservative quantities of the evolutionary equation. To 
distinguish this class of symmetries from the previous one we will call them scaling symmetries. We show 
that the restriction of the starting equation on the finite dimensional manifold admits a natural description 
as a Hamiltonian system with time-dependent Hamiltonian. 

The best known class of examples of evolutionary PDEs admitting nontrivial symmetries and 
conservation laws are integrable systems of soliton theory (see [SM] and references therein). The finite 
dimensional manifolds of the stationary points of integrable systems are typically described by ODEs of 
Painleve type [AS],[CD]. For the simplest examples of these restrictions the Hamiltonian structure is already 
known. For example for the classical six Painleve equations the Hamiltonian description was found by 
Okamoto, [O]. Although the relationship between the starting PDE and the reduced ODE is clear and has 
been investigated quite a lot (see, e.g. [AS]), the relationship between the starting Hamiltonian structure 
and the reduced one has not been elucidated. This work will give a contribution in understanding of this 
relationship. 

As a first result (see section 2) we prove that the finite dimensional Hamiltonian structure of the 
ODEs is obtained from the Hamiltonian structure of the starting PDE, via scaling reduction. Particularly, 
we construct the time-dependent Hamiltonian function of the reduced system. In the time-independent 
case this procedure coincides with the well known stationary-flow reduction discovered by Bogoyavlenskii 
and Novikov [BN]. As an application we present the case of PI, PH, PHI, PVI and also certain higher order 
systems appeared recently in the theory of Frobenius manifolds [Dl]. 

As a second result (see section 3) we present a very general Lagrangian formulation of the procedure 
of reduction of an evolutionary system (0.1). Namely, we prove that this restriction is again a Lagrangian 
system with the Lagrangian function A, such that 

dk dL 
dx dt 


The work is structured as follows: after recalling, in Section 1, some basic facts about the Hamiltonian 
structure of the evolutionary PDEs, and briefly summarizing the method of reduction of evolutionary flows 
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on the manifold of stationary points of their integral, introduced by Bogoyavlenskii and Novikov [BN], in 
Section 2 we consider the generalization of this procedure to scaling symmetries. The reduced flow is a time- 
dependent Hamiltonian system, and in Theorem 2.1 we give the relationship between the infinite-dimensional 
Hamiltonian structure and the reduced one. 

Section 3 is devoted to a Lagrangian approach to the problem: after describing the general framework, 
in Theorem 3.1 we give the procedure of reduction and we construct the reduced Lagrangian function. In 
Section 3.2 we establish the relationship with the Hamiltonian approach. As an application we study the 
Lagrangian reduction of KdV on the manifold of the fixed points of the 7-th flow. 

Section 4 contains the application of the theory to the scaling reductions from KdV, mKdV and Sine- 
Gordon equations respectively to Painleve I, Painleve H and HI. These examples are studied both from the 
Hamiltonian and the Lagrangian point of view. 

In Section 5 we study the n-waves equation and his scaling reduction to a system of commuting 
Hamiltonian flows on the Lie algebra so(n). The reduced system is a non-autonomous Hamiltonian system 
w.r.t. the Poisson structure of so(n). In particular, for n = 3, adding an additional symmetry condition, 
one arrives at Painleve VI equation. 


1. Infinite dimensional Hamiltonian strnctnres and stationary fiows reduction 

Let us consider the phase space dJl of smooth maps of the circle into some smooth n-dimensional manifold. 
Actually we can forget about the boundary conditions when dealing with local functionals only. We denote 
by U the space of smooth functionals on Tl of the form 

F{u) = J f{x,u{x),Ux{x),... ,u^'^\x))dx, 

where the density / depends only on a finite number of derivatives of u. On the space U the variational 
derivative is defined by 

r 6F 

SF = / 6u^(x)dx. 

J (5fo(x) ^ ’ 

Explicitly, 

^F ^ df ^ y. / df 

6u^{x) du^{x) ^ dx'" du'F)[x) 

One can define on 911 the (formal) Poisson brackets 


N 


{u\x), (y)} = w' (x, y)=Y^ Afc d'"’ {x - y) 




where Al depends on a finite number of derivatives of u. This induces on ^ the Poisson bracket 


{F,G} = 


6F ii 6G , 
--r , ,. , ax, 


Su'^{x) Su^{x) 


where 


N 






A Hamiltonian system on DJI has then the form 


l{x)={u\x),H} = P^ 


SH 

5u^ (x)' 
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In particular, we consider so called Gardner-Zakharov-Faddev bracket P' 
Hamiltonian system has the form 

utix) = {u{x),H} = ^ 


ax 


with Poisson bracket 


Let us consider a first integral 


dx 5u{x) 


(F,G}= 


( 1 - 1 ) 


Su{x) dx Su(x) 


/ = y L{x,u{x),Ux{x),... ,u^'^^)dx, 
where L does not depend on t. The generalized Euler-Lagrange equation 

SI 


Su{a 


= 0 


( 1 . 2 ) 


generically is a ODE of order 2n fixing the 2n-dimensional manifold 6 of the stationary points of the first 
integral I. Because of the Lax lemma (see [Mo]) this submanifold is invariant w.r.t the evolutionary equation 
(1.1). The functional L is the Lagrangian of the x-flow defined by (1.2). If L is nondegenerate, then it defines 
also on 6 the natural system of canonical coordinates 

t = l,2,....,n 

- JL 

and equation (1.2) can be put in the Hamiltonian form 


iP.h = 

= Wp 


where H is the generalized Legendre transform of L\ 


h = -l + Y. 


51 

6u^^'> 


,(b 


which, in terms of the canonical coordinates takes the form: 


h=-l+y: 

1 



It is well known that the starting PDE can be restricted on © and the restriction is a Hamiltonian system 
of ODEs. In particular Bogoyavlenskii and Novikov discovered the algorithm to construct the Hamiltonian 
functions of the reductions in terms of the Hamiltonian of the original evolutionary equation. They considered 
the case of a hierarchy of evolutionary equations 

du d 51,^ 

dt^. dx 5u{x) 

with = J {u{x),Ux{x ),..., {x))dx, and they described the reduction procedure of the fc-th flow on 

the finite dimensional manifold of the stationary points of the j-th flow. They proved that all the flows of 
the hierarchy reduce to finite dimensional Hamiltonian system. The Hamiltonian function for the reduced 
/c-th flow, {—Qk,j), is determined by : 
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d SI, 

6u{x) dx Su{x) dx 

Mokhov [Mo] generalized this result to not necessarily Hamiltonian evolutionary PDEs. 


2. Scaling reductions of evolutionary systems: 
Hamiltonian formulation 


In this Section we extend the Bogojavlenskii-Novikov scheme to finite dimensional invariant submanifolds 
specified by time-dependent local symmetries. 

We start from a partial differential equation of order m on the functional space describing the 
evolution of the function u{x) in the time t and a scaling symmetry 

Ua = G{x, t, U,Ux,..., 

Our main assumption is that the set of stationary points of the symmetry can be formally represented in 
the Euler-Lagrange form 

5u{x) 


/ = y L{x,t,u,Ux, ■ ■ ■ ,u^"^'^)dx, 


It is an ordinary differential equation of order 2n depending explicitly on the parameter t. If L is 
nondegenerate, the space of the solutions is a 2n dimensional manifold ©, which naturally carries a system of 
canonical coordinates. As in Section 1 we will show that, in these coordinates, the Euler-Lagrange equation 
is Hamiltonian, with Hamiltonian function H, obtained from L via Legendre transform: 

TT r I 

i=l 

Following the scheme of [BN], we prove that one can reduce on © also the equation of the evolution in t, 
which results to be a Hamiltonian system. We also give a universal scheme to produce the time-dependent 
Hamiltonian function of this reduced system. Indeed the following theorem holds: 


Theorem 2.1: If the evolutionary PDE: 


= F{u,Ux ,.. 


admits a nondegenerate scaling symmetry, then, on the manifold © of the stationary points of the symmetry: 


1= L{x,t,u,Ux, ■ ■ ■ ,u )dx, 


it reduces to a Hamiltonian motion in t, for the time dependent Hamiltonian function {—Q), that is the 
reduction on © of 

Q = (2.1) 
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( 2 n — 1 ) 

where Pi, qi are the canonical coordinates on &, expressed in terms of u, u^, ■ ■ ■ ,u 
is determined by 

dL dA 
dt dx 


and the function A 


( 2 . 2 ) 


Proof: We prove the theorem in three steps: first we describe the submanifold 6 of stationary points of 
the symmetry I, where we introduce a system of canonical coordinates; then we deduce, on ©, a zero- 
curvature equation for {—Q) and the Hamiltonian function H of the reduced a;-flow. Finally we prove that 
the restricted t-flow is Hamiltonian on ©,with Hamiltonian function (—Q). 


1) The manifold © is the 2n-dimensional manifold of the solutions of the Euler-Lagrange equation 


SI 

Su(x) 


= 0 . 


(2.3) 


It is invariant under the foflow and it naturally carries a system ofcanonical coordinates: 


g. = ttb z = 1, 2,..., n 

- JL 


(2.4a) 

(2.46) 


obtained via generalized Lagrange transform (here we suppose that the generalized Lagrangian L is 
nondegenerate). Observe that now the pi depend on x and on t. 

Reversing relations (2.4), one can express the derivatives u,Ux, ■ ■ ■ in terms of the canonical 

coordinates pi and qt, x and t; explicitely: 


= {qjx = giix,t,q,,...,qn,Pn) 
■uO+l) = g 2 {x,t,q^,...,qn,Pn,Pn-l) 


■u(2" 1) = g„{x,t,q^,...,qn,Pn,---,Pi)- 


Observe that (2.4) gives the identities: 


C„ 'I _L cm ^ 
yPlJx -I- — Su 

(Pi)^ + if- = 0, * > 1 


(2.5) 


(9.).-|f ^0, 

where H is the generalized Legendre transform of L: 




SI 

Su^'^') 


,b) 


which, in terms of the canonical coordinates takes the form: 




dx 


The first of identities (2.5) allows us to express the higher derivatives for m > 2n in terms of x,t, 
Pi, qi and pi^ with ? = 1,...,m — 2n + 1, explicitly: 

= gn+l{x,t,q^, . . . ,qn,Pn, ■ ■ ■ ,Pi, {Pl)x) 


_ / , / \(^ 2 n + l) , 

^ — 9m—Qi Qri'! Pn-) t Pii 1 [Pi ) )• 
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On 6 it reduces to {Pi)x + ^- = 0, and the system (2.5) is a canonical Hamiltonian system, with Hamiltonian 
function H, giving the reduced a;-flow. 

Now we will show that also the t-flow reduces on 6 with Hamiltonian function (—Q). 


Firstly we observe that Q is a function of x,t, u and its x-derivatives up to the order (m + n), then 
it can be rewritten in terms of x, t, (pi, qi) and p® up to the order I = m — n -\- 1. 

We denote with / a function f{x,t,u{x),... ,u{x)^^'>) reduced on 6; notice that, if j > 2n, then the 
reduction can be done using the relation 


iPl)x 


dl^ 

dq. 


( 2 . 6 ) 


Then / does depend explicitly only on the pi and qi, for i = 1,... ,n and on the time t. In fact differentiating 
(t.6) one obtains the derivatives p® in terms of the canonical coordinates {pi,qi). 


2) We consider the derivative 


dL dL ^ dL dqt ^ dL dpi 
dt dt ^ dqi dt ^ dpi dt 

l — L l — L 

^ d?q, 

~ ~ ~ ^ 7)^ TiT ^ 


dqi dt ' ^^‘dxdt 
1—1 1—1 


(2.7) 


From the fact that / is a first integral, one deduces that ^ must be the total derivative in x of a functional 
A that does depend on x,t,{pi, qi) and pi'’ up to the order ^ = m — n + 1; we have: 


dA _dQ ^dQdq, 


\" dQ dpi f \{i) 


dx dx ^ dqi dx ^ dpi dx ^ dip, )d) dx 

i—1 i—1 i—1 -1- / 

dpi dqi d7qi _ 

^ dx dt dxdt 

2=1 2=1 

dQ d (,) 

r)rr / j firi: fin- f ftrt- fin- ‘ \(i\ J™ 7^1 


^ (9(pjW<ia; 

dpi 


dH dqi dqi . , dQ . 

“ + —(Pi)- 


dqi dt ' ^^"dxdt ' dt 

i=2 ^ i=l 


Then equation (2.2) gives: 


dH ^ dQ 

dt dx dqi dpi ^ dpi dqi dpi 


m—n+l , , 

y- dQ d dqi ^ 

^ d{pQd)dx^^^’ ^ dt^^^' 


dH dqi 
dqi dt 


= 0 , 


which can be rewritten as 


dH dqi 


dt dt 


dH 


+ ^ ipQ^ + 1^ =-^Q 


dqi 


dx 


At this point we need the 

Lemma 2.1: On the submanifold & the following relation holds: 


dipQ<~J) 


= 0 Vj > 1. 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 
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Proof: See Appendix 2.A 


Hence, on the submanifold 6 eq. (2.9) reduces to: 


dH dQ 
dt dx ' ^ dqi dpi 


^dQ dH ^ dQ dH 


i=l 


dpi dqi 


This is a zero-curvature equation: 


This completes the second step in the proof of the theorem. 

3. Now we will construct the Hamiltonian system inductively; to this end we need a further lemma: 
Lemma 2.2: The fundamental relation 


dqi 

dt 


dQ 

dpi ■ 


( 2 . 11 ) 


holds. 

Proof: See Appendix 2.A 

For simplicity, here and in the following we omit the “tilde” sign: Q will indicate the reduced function 
on ©. Now, we assume that ^ Q} and we prove inductively that the same relation holds 

for q^_^_l. The scheme of the procedure is the same as in [BN], the only differences are the contributions of 
the partial derivatives in t and x. Indeed, 

Using the Jacobi identity and the zero-curvature relation we get 

dQi+i 




z = 1, 2,..., n — 1. 


Here the term {^,g.} is zero for every i ^ n since (^) = — Indeed L depends on u and on the 
derivatives of u up to the order n. This means that, restricted on ©, it depends on q^,q^,..., qn+i■ Then, 
there is no dependence on the p^ for i ^ n. Finally we get 

dH 

j Qn} T Oj 




I < n. 


Hence we have proved that 


d(li 

dt 




i = 1,2,n. 


Now we prove that = 


— {Pi, Q} by induction, starting from p„. 



This comes from the commutativity of the flows 


d d 


d d 


dt dx^'^ dx 


explicitly: 


-(—( 

dt \dx 


Qn I — 


__d _ 

dt V dpn ) 

d dH ^ dH dQ d‘^H d 


i=l 


dpndq^ dp^ dpl^ dr 


dH 


=fe.^)-E 


dH dQ d^H d 


7Pn- 


dt ’ ^ dpndq^ dp^ dp^ dt 

On the other hand, using the Jacobi identity and the zero-curvature equation, one can write 

d dQ 


d f dq„ 
dx \ dt 


dx dt 


- {{qn,Q},H} = 


do 

= - Un, {Q, H}} - {{qu,H}Q} - {qn, -^} = 
dH dQ d‘^H dQ 






dpndq^ dp. dpi dqr, 


Comparing the two expressions and noticing that 0 because of the nondegeneracy, we get 

dpn r dQ 

IF = -<""■« = 

- Now we suppose that ^ = {p., Q} and we deduce the same for p._^. Indeed, 


d / dp^ 
dx \ dt 


^ dxdq, 

= -{{p^,Q},H]-{p„^} = 

f) M 


and 


=fa.^)+i: 




dt \ dx 


dt \ dq^ 

d dH 

dt, dq, 


d dH ^ 

f)-f r)n < ^ 




i=i 


-4 d'^H dQ 

i - 

^ dq^dq^ dp^ 

d‘^H dQ 

dq.dpn dqn 

d'^H dQ 

dq^dp^_,i dq^_^ 

^ d'^H dq^ 

d‘^H dpn 

d‘^H dp._.^ 

^ dq,dq^ dt. 

dq^dpn dt. 

dqidp^_^ dt^ 

d^H dQ 

d'^H dQ 

d^H dp^_.^ 

dq^dq^ dp. 

dq.dpn dqn 

dq^dp^_.^ dt, ’ 
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where 


d'^H 


1. Comparing the two expressions we get 


dPi_i _ dQ 
dt 9Qi_i ’ 


hence it follows that 


dp. dQ 
dt dq^ 




i = 1,2,..., n 


Q.E.D. 


Remark: The definition of Q: 


-Q 


-a + E 



looks very similar to the definition of the Hamiltonian function H of the x-flow: 


h = -l + Y. 

i=l 



Here a symmetry between x and t seems to appear: one could be tempted to read the definition of Q as 
a Legendre transform and hence to read A as the Lagrangian of the t-flow. But it is not completely true: 
indeed the coordinates qi and pi are obtained from the Lagrangian L, they are not, a priori, good coordinates 
for A. In the next chapter we will perform a change of coordinates on ©, in order to read A as Lagrangian 
function. 


2.A Appendix 

Proof of Lemma 2.1: We observe that the recursive relation 

{ dQ \ d f d \ d / dQ \ 
\d{pQU-^)) d{pQ(d \dx^) dx\d{pQ(d ) 

holds for j > 1. Indeed 


d 


d{pQ^d ydx 


m—n+1 

+ E 


d^Q 


dq^d{pQ(d 

d^Q 


{di)x + 






ipt)a 


2=1 


d{pi)^dd(^pQU) 




dQ 


d{pQU-^)J- 


When we reduce on ©: 

f dQ \ ^ 

\dipQU-^)J 

But g depends on {pQ^^^ up to a finite order, then 


d / dQ \ 
dx \d{pQ^d J 


dQ 


d{pQ(d 


0 Vj > 1. 


Q.E.D. 


Proof of Lemma 2.2: The expansion of -^Q in powers of {pQx near the point reads 



~ d 

APi) 



(Pl )x + 


dH 

dqi 


0 (Pi)a 


dH 

dqi 
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The zero order term is 




dx 


dt 


by virtue of the zero-curvature equation. The first order coefficient is 


dipjx \dx 


-Q 


d { d 


dx ydipj: 

d f ^ 


dx 


+ 




m —n+1 

k 


dQ 

dpi 


where the only non zero term is by virtue of Lemma 2.1. Hence we obtain the power series 
expansion of -^Q up to the first order: 


dH dQ 


dt dp^ 


dH 


-^ +IT- (Pi)- +TTT 


dqi 


this, compared with the left-hand side of eq. (2.9), gives the relation (2.11). 


Q.E.D. 


3. Scaling reductions of evolutionary systems: 
Lagrangian formulation 


3.1 General framework 

The basic idea is to develop a reduction method dealing on the same footing with x and t. The starting 
point is always the evolutionary PDE 


Ut=F{u,Ux, .(3.1) 

in the space fDl described in Section 1.1. The first step of our construction is to read u as a function of x 

and t and to consider equation (3.1) as a definition of in terms of u{x,t),Ux{x,t), . ,u^"^~^^x,t) 

and Ut(x, t). 

This corresponds to consider as “coordinates” in 911, instead of u{x,t) and its derivatives in x: 

U, Ifa;, ‘Uxx'i ■ ■ ■ 

(here and in the following u indicate the function u{x,t)), the system 

('777_i'\ (m—1) 

U,Ux,...,F \Ut,Uxt,---,ul 

By virtue of the reversibility of (3.1) in u^'^\x,t) it is possible to perform this “change of variables”. If one 
introduce the vector 

u = {u,Ux,. ■. 

the new system of “coordinates” in dJl is given by u{x, t) and its derivatives in t: 


tt, Uf) ^tt'i • • • 
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At this point one takes a first integral of eq (3.1), i.e. a functional 


I = 


L[ x,t,u{x),Ua;ix), 



dx 


(3.2) 


in the space Tl, such that 


61 

Su{x) 


= 0 . 


(3.3) 


This Euler-Lagrange equation defines a finite dimensional manifold ©, i.e. the set of the fixed points of 
I. Indeed the Euler-Lagrange equation (3.3) is an ODE of order 2n, so that the space of the solutions is 
a 2n-dimensional manifold; © is modeled on this space, having as coordinates certain combinations of the 
initial values, i.e. of the first (2n — 1) a;-derivatives of u{x) evaluated at xq. 


In Lemma 3.2 below, we rewrite the definition of the manifold © in terms of u{t), ... and of a functional 


J = 


A^x,t,u{t),ut{t),... 


where A can be calculated from L (see eq. (3.4)), and the order [3 of derivation in t depends on the ratio 
between m and n, as we will show in detail in Section 3.3. 

In Theorem 3.1 we will prove that A^a;, t, u{t), Ut{t), ..., is the generalized Lagrangian for the t-flow 

reduced on ©. Indeed equation (3.1) can be rewritten in form of a Euler-Lagrange equation: 


dJ 

Su{t) 


= 0 , 


for the vector 


where 


6J / 6J 6J 6J \ 

Su{t) \Su{t)^ Sux{t)^ ’ 


SJ 


dA 

i9u(*.o) 




dA 


In the multiindex (t, a) the Latin character indicates the order in the x-derivative, the Greek indicate the 
order in the t-derivative. 

Explicitly, equation (3.1) reads 


dA 

dA 

dux 


E(-i)“ 


dA 

dt^ du^°'') 

dA 


0 

0 


dA 


E(-i)“ 


d“ dA 
dt" 


0 . 


We will formalize these facts in the following Theorem 3.1; here we give an idea of the proof, ignoring all the 
calculations, that we will concentrate in Lemma 3.1 and 3.2, the proof of which is postponed in Appendix 

3.A. 

The proof is performed in four steps: firstly we define the new system of “coordinates” in 911, and 
we give some useful relations between the new and the old “coordinates”. As a second step we rewrite the 
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Lagrangian density L(u{x)) in terms of the new “coordinates” and we construct, starting from L(u{t)^, the 
Lagrangian density A(u{t)). 

The third step consists in recovering the relation between performing the variation of L (in x) and of 
A (in t). The most relevant relation is that given in Lemma 3.1. This relation is necessary to rewrite the 
Euler-Lagrange equation (3.2), defining ©, as a condition on A. The explicit form of this condition is given 
in Lemma 3.2. 

Finally we prove that, under this condition, i.e. after performing the reduction on ©, the starting 
evolution equation (3.1), reads as an Euler-Lagrange equation for A. 

The method of Hamiltonian reduction described in Chapter 2 allows us to put a canonical system of 
coordinates {pi, qi} on © (see formula (2.4)). These coordinates are obtained from L via generalized Lagrange 
transform, so that they are, in a certain sense, adapted to the x-flow. This means that in these coordinates 
the reduced x-flow is a Hamiltonian system. Theorem 2.1 also gives the explicit form of the Hamiltonian 
function 

H = -L +y^Pi(gi)x- 

i 

The method of Lagrangian reduction which we describe in this Chapter, still allows us to define a system of 
canonical coordinates: we will call it {pi,qi}. These coordinates are obtained from A, i.e. they are adapted 
to the t-flow; in fact we will prove (see Section 3.3) that, in these coordinates, the reduced t-flow is a 
Hamiltonian system, with Hamiltonian function 

-Q = -A + '^p,{q,)t. 


When rewritten in terms of {pi,qi}, the Hamiltonian Q coincides with the Hamiltonian function Q 
constructed by Bogoyavlenskii and Novikov. 

In this sense the alternative definition of Q given by us in Theorem 2.1: 

-Q = -A + '^p^{q^)t, 

i 

is a Legendre transformation, if one uses the right system of canonical coordinates (see below). 


3.2 Lagrangian reduction 

Theorem 3.1: If the evolutionary PDE: 

Ut = F{u,Ux, .(3.3) 

admits a nondegenerate scaling symmetry, then, on the manifold © of the stationary points of the symmetry: 

61 


Su{x) 


= 0 , 


/ ( ) dl 

L{x,t,u,Ux,...,u )dx, 

it reduces to a Lagrangian motion in t, for the time dependent Lagrangian function A, determined by: 

dL _ dA 
dt dx 


(3.4) 


13 




Proof: (in the following we consider the case m < n < 2m. The same holds in the case (a — l)m < n < am, 
as we will show in Section 3.3). We prove the theorem in four steps: 


1. Change of “coordinates”: Let us assume that the evolutionary equation (3.3) depends on u(x) and on 
its x-derivatives up to finite order m, and that this equation is invertible in In this case we can read 

(3.3) as a definition of in terms of u, Ux,..., and Ut- 




(3.5) 


Differentiating eq. (3.5) in x one obtains all the x-derivatives of u of order greater then m in terms 
of u, Ux, ■. ■, and their t-derivatives: 

f = fi{u,Ux, ■ ■ ■ ,U^"'~'^\ut,Uxt) 


= fm{u,Ux,...,U^"^ '^\ut,Uxt---,Vf ,Utt) 


,Ut,Uxt ■ ■ ■ ,Ul ') 
m — 1) 
t 


[ = /„(m,'U2:, . . . ^\ut,Uxt---,u\"' ^\uu,Uxtt, ■ ■ ■ ,U^u ’"^)- 

Explicitly, the first relation has the form 


(m+l) _ 


+ 


-Ux +...+ 




,(™) 




dx ' du ' ' du ^”“-~'^">" ' dut 

and in general, using the multiindex notation introduced in Section 3.1, 






dx 


m—1 a 

EE 


5y(m+j-l) 


,(fe+l./3) 


fc=0 / 3=0 

where the higher order a in the t-derivative is fixed by {a — l)m < j < am. 


This completes the construction of the map from 


U,Ux,..., . . . , . . . , „(2m)^ 


to the new system of “coordinates” 


(m — 1) (™ —1) (n-m) 

U,Ux,...,C >,Ut,...,Ul \Utt,Uxtt,---,ult %. 


Here below we list some noteworthy relationships between the two system (they will be useful in the 
following): 

(„) _ du^ duP^ („_i) 


u; = 




duf 


-Utt, 


(3.6a) 


d /agWX 
dx \ dut J 
d / aaW \ 
dx\dui’'W 
d / aa(^) \ 

dx \ du^C j 
d ( du^^'> \ 

dt \ dut J 


au(*+i) au(*) au(™) 

dut dut 

au(*+i) aw^*) 

au(®+i) aw^*) du^‘^ du^"^^ 

du^C du^^~^'> du^C 

a^ _ a^ 

dut du 


(3.66) 

(3.6c) 

(3.6(i) 

(3.6e) 
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2. Lagrangian densities: The Lagrangian L defining the symmetry, depends on u,Ux, ■ ■ ■ so that its 

derivative ^ depends on u,Ux, ■ ■ ■ In terms of the new “coordinates” one may rewrite L as 

L{x, t, u,Ux,..., ut, ■ • ■, (3.7) 


and 


dL\ fdL\ fdL\ 

=hi 7 +h;-Wt + --- + 


dt J \dt J \du J 


dL 


J t 

'dL 


(m—l) 




Of course, (3.8a) coincides with 


dL dL dL 


aL 


au(m-l) 


(m — l) 


aw, 


dL 


■Utt 




(n—r: 


where 


^-(^\ V 


k—m 


dL 


dL 


duf'd V du^d 

dL 


\du(>^)) dt 


^ I aitW ] au(a 

k—m 


du. 


id 


= E 


/ dL \ aaW 

Vau(o j goa 


t = 0 ,..., m — 1 

z = 0,..., TO — 1. 


From the fact that I is a first integral, it follows that there exists a functional 

A Z J- (m—l) (™ —1) (n—m —1 )n 


such that 


where 


dL _ dA 
dt dx ’ 


dA _ dA dA 
dx dx du 
dA 


dA 


a„(m-2) 




[U,Ux,...,U ,Utt)+ g^Uxt + ■■■+ 

dA (jji'j , , dA dA (n—m) 

+ {u,..., Utt) + ^Uxtt + . . . + «ZZ ■ 


du\. 


(3.8a). 


(3.86) 


(3.9a) 

(3.96) 

(3.9c) 


(3.10) 


3. Variations: Our aim is to reduce equation (3.3) on the space © of the stationary points of I = f Ldx. 
This finite-dimensional manifold is defined by the Euler-Lagrange equation 


5u{x) 
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This is a variational equation in the old “coordinates” u,Ux, ■ ■ ■] how can we define the same manifold © in 
terms of the new “coordinates”? We must express in terms of A and its variations. To this end we 
recall that 


SI 

6u{x) 


— + Vf-lV — 

du ^ dP du^P' 


and we first express the terms 


5L 

Su^P{x) 


for j > 0 in terms of A and the the new “coordinates”, namely: 


Lemma 3.1: The following recurrence relation holds: 


dk d f dk \ 
V du[]^ ) 



/ SI \ 

J duPl 


i = 1,..., TO — 1 


(3.11) 


Proof: see Appendix 3.A 


The proof of Lemma 3.1 is based on the comparison of (3.8) and (3.10) and their partial derivatives w.r.t. 
and u[{\ With a similar technique, and using equation (3.11), one can proves the fundamental 

Lemma 3.2: The (generalized) Euler-Lagrange equation 


Su{x) 

is equivalent to the condition: 

dk d f dk 


= 0 . 


(3.12) 


Proof: see Appendix 3.A 
Introducing the functional 

J = y A ^a;, t, u{t), Ux{t),..., (t), ut{t),..., (t)^ dt, 


equation (3.12) reads 

=0 

Su(^-P{t) 

Notice that the object in the left hand side is the last component of the vector 


SJ _ / SJ SJ SJ \ 

Su{t) \i5u(t) ’ Sux{t) ’ ’ ’ Su^'^~P{t) J 


4. Reduced evolutionary equation: Here we prove that all the components of the vector are zero. This 
can be done recursively, by mean of 
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Lemma 3.3: The following recurrence relation holds: 

SJ 6J d f 5J \ 

(5u(*-i)(t) dx\6u^^l{t)) * 


(3.13a) 


Proof: see Appendix 3.A 

Indeed, Lemma 3.2 states that the (to — l)-th component is zero when reduced on ©, hence, by 

virtue of (3.13a), all the components of vanish on 6. 

This is the Euler-Lagrange equation for the Lagrangian 


A t j. (”“- 1 ) 

A(x,t,u,Ux, ■ ■■ ,u 




(m-1) 


Utt, ■■■,U 


(n—m) 

tt 


Q.E.D. 


Remark: Equation (3.13a) can be rewritten as 


■ dA 


f dA \ 

\+^i 

f dA Y 


dt\ 



Kdu[l~^U. 


dA 




dA 


1) dt\dJ'^-dJ\ 


dul-d 


A. 

dx 


■ dA 


^ dA " 

U 1 

f dA Y 

du^d 

dt ' 





for i = 1,...,TO — 1. And The Euler-Lagrange equation reads 


(3.135) 


3A 

duO) 


d dA _ r. 


i = n — m, 1 


dA _ d ( dA \ I d^ dA 
du'-') dt V ' dt'^ 


duf- 


TT 


= 0 


I = 0,..., n — TO — 1. 


(3.14) 


3.3 Relation with the Hamiltonian reduction 


Theorem 3.1 provides an alternative definition of the space ©, and of the relative system of canonical 
coordinates: 

i = 1,..., TO 


q,=q,= 


'Im+i = (9i)i = 

Pi = 


^ i = 1,... ,n — m 

jtAtt- z=1,...,to 


Pm + i = 


dA 


(3.15) 


du 




i = 1,..., n — TO 


We consider now the Hamiltonian (—Q) of the reduced t-flow, defined in Theorem 2.1 


Q = A - y^Pi{qt)t = A - 


n—1 


SI 


,(i) 


i=l 


^ (5ub+i)(t) 

2=0 ^ ^ 


At the end of the Chapter 2 we noticed how this expression looks very similar to a Legendre transform, but 
it is not; here we will show that actually the Legendre transform of the Lagrangian A gives the Hamiltonian 
Q, where Q is written in the coordinate system relative to A. 
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Firstly we rewrite Q in the coordinate system (3.15): 



Using Lemma 3.1 we get 


m —1 c T n—m—1 o* 

A * X —' (i) \ —' uA. 

-O'-A+ETTTyd” 

i=0 H 


E 


o (d “ 
i=o 




for A(x,t,qi,{qi)t,...,{qn-m)t)- 


(3.16) 


3.4 Concluding remarks 


The case considered in Theorem 3.1 is the more general one. Indeed, for {a — l)m < n < am, the Euler- 
Lagrange equation defining ©: 

'61 


Su{a 


= 0 


into the new “coordinates”, is a differential equation in u,..., • ■ •, . 

The Lagrangian L transforms into 

L{x, ..., uf^~^\uu, • ■ •, 

and we can define the new Lagrangian 

A/ I (m—1) (^“1) (n—'m—l,a)\ 

A(x, t, U, Uxj • • • ,U , Ut, . . . , , Uttj • • • jU ). 

The proof of the theorem is the same, one has only to consider the identities 


d fdL\ _ d fdA 


(9u(*./3) \dt ) du^'^’d) I dx 


for i = 1,...,TO — 1 and P = 1,... ,a. 

In particular, the manifold © is defined by 

=0 

(5'u(™-i)(d 

and it naturally carries the canonical system of coordinates 

{ '?/ 3 m+i = I = 1 ,...,to; P = 0,... ,a — 2 

Q(c,-i)m+i — i = 1,... ,n — {a — 1)to 


(3.17) 


(3.18) 
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In the following we will consider in detail the case 0=1, i.e. n < m, which occurs in the applications we 
are interested in (see next chapter). In this case L and L coincide and 


dL dL dL dL („) 

-ut + ... + ■ 


dt dt du 




The new Lagrangian is 




with 


dL _dA _ dA dA 
dt dx dx du 


Ux + ■ ■ ■ + 


dA f ) dA dA („) 




du 


(n-1) 


(3.19) 


(3.20) 


In this case Lemma 3.1 reduces to the following recurrence relation: 

SI dA 


6u^d(x) dul" 

and the proof is based on the identity 


i = 1,...,n 


d f dL\ d f dA 


dui"^ \dtj V dx 


— i = 0,...,n-l, 


observing that, for i>l, 


d fdL\ _ dL 


Qy^) \dt J du(d ■ 


In fact, L does not depend on u[^\ and 


d d ( dA 


dA 


(3.21) 


(3.22) 


duf^ \dx) dx\ auW ) duf ’ 

In particular the first step, i = n, follow directly from the fact that the only dependence of in both 
(3.19) and (3.20) is the one explicitly shown, so that 


dL 


dA 


On the other hand, from (3.22) for the index z = 0, one obtains the fundamental relation 

61 _ dA au(™) 

5u{x) dut 

and, since is always nonzero, the condition that defines the submanifold © is 

au(m-l) 

The relative system of canonical coordinates is given by: 

S n = 9A 

1 


(3.23) 


(3.24) 


For i = 1,..., n. 
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The reduced t-flow is Lagrangian, with Lagrangian A. 
Indeed, from the identity 


^ ( 

UL\ 

d 


du^^) ' 

V dt ) 

i9Mb) 

\dx J 


i = 1,..., TO — 1, 


and using the Lemma, one obtains on the subspace ©: 


aA 


= 0 


a«co " * n,..., TO 1 


* = 0,...,n-l 


which is the Euler-Lagrange equation for the Lagrangian A(a:, t,u,Ux, ■ ■ ■ ,u 


(m-l) 


,Ut, 




3.5 Example: KdV with tr fixed 

We will give below an example of how does Theorem 3.1 works for the first non trivial case, n = m. We 
study the Lagrangian reduction of the KdV equation 

ut = 6uux - Uxxx (3.25) 

on the stationary manifold of the ty-flow. 

The Lagrangian density of the t 7 -flow, reduced to the normal form (here I mean that L does not 
contains total derivatives), depends on the a;-derivatives of u{x, t) up to order n = 3, and has the expression 

L = 7u^ + 35u^ul + 7uulx + ^(u^^^)^ (3.26) 

The submanifold © of the stationary points defined by the Euler-Lagrange equation for L gives the n+m = 6 
derivative in terms of the first five, explicitly 

.j^(6) _ i4yy(4) _|_ 2SuxUxxx — TOu^Uxx + 21w^a; — 70uul + 35m‘*. (3.27) 

From the relation 

dL dA 
dt dx 

one construct the Lagrangian A(u, Ux, ■ ■ ■, By direct calculation 

A = - f - lAuUxxU^"'^^ + + l^UxUxxUxxx + 

—70u^UxU^^^ — (m^^^)^ + 77u^{uxxY + 70uUxUxx — 35u‘^Uxx+ 

35 

+ 280u^ul + 35m®. (3.28) 

The evolution equation (3.27) is the definition of Uxxx in terms of {u,Ux,Uxx,Ut), explicitly: 

^xxx — 6MMa, Ut- 


Differentiating this relation in x one obtains 

{ ^( 4 ) = 6uuxx + 6ul - Uxt 

m(5) = ISUxUxx + 36m^ - 6MMt - Uxxt 

^(6) = ISm^,,, + 180mm^ + 36u^Uxx — 30uxUt — 12uuxt + Uu 
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which is a map from the “coordinates” 




U, llx-i Uxxj '^xti '^xxt-, • ■ ■ 

The Lagrangian L depends on u,Ux,Uxx,Uxxx] in the new “coordinates” 

L{u, Ux, Uxx, ut) = 7u^ + + Tuu^^ - duUxUt + 

Its derivative looks like 

at 

dL di di di di dL 

“77 — “TTT I a "o '^xt H" 7^ ^xxt H" "7 '^tt- 

dt ut OU UUx UUxx UUt 

And there exist a functional A depending on u^Ux^Uxx^ut^Uxt^Uxxt^ explicitly 
A =35^® + iu^Ux + -Ux — dUxUxt + 

35w Uxx ‘^UUxUxx “t“ ^UUxfUxx “t“ llw Uxx~^ 

-ulx + GuUxUxxt + 22u^UxUt + AUxUxxUt - UxxtUt + ^UUt 


such that 


where 


dL dA dA dA dA dA /o-, 

—JT ~ T Ux H Uxx “t“ "T U 

dt dx ax OU aux auxx 


dA dA 

Uxt “t“ 

dut duxt 


dL A o o 

= 35u^ + 106 mu^ + 7u‘i„ - QuxUt 
du 

—— = 106u‘^Ux — 6uut 


dA 


- — ^'^UUxx 

jUxx 

dL 

-— = -6uux + Ut 
UUt 


= 2l{)vr + 12u^u^ — lAQvruxx — 2u^Uxx + SuxtUxx + + ^UxUxxt + 4:AuUxUt + 

= Su Ux 2u^ \2%lxUxt ‘^UUxUxx ^UUxxt 22u Ut ~\~ AiUxxUt 


= — 35 ^^ — 2uul. + ^uuxt + 22u^Uxx — ^u^^ + ^UxUt 


— 22u Ux ^UxUxx Uxxt “t“ Suut 


= -6ui + Uxt + Suux 


=6uux — Ut 
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Lemma 3.2 states that the condition 4^ = 0, which defines the submanifold ©, is equivalent to the condition 


dA _d_l' dA 


du 

XX dt Y diixxt 


= 0 , 


explicitly: 


Utt = 35u^ + 2 wm^ — 22u^Uxx + ‘^u’tx + ‘^UxUt — 2uuxt- 


The Euler-Lagrange equation for A reads 


, (^ 35w 2uUx T 22u Uxx ^^xx ‘^'^x'^t T 2,uuxt — 0 


du 

XX dt diixxt 

dA d , dA . o „ , r. ^^2 

— ”77 ("o ) — ^^x ^^^x^xx ^'^'^xxt 221/. Ilf Allxx^t ^xtt — 0 

dux dt ^ duxt 

dd / d \ c: O O Q O O 

-17(-^ —) = 210u^ + 12ia - lAOu^Uxx - ‘^u^Uxx + ^UxtUxx + 22itit^^+ 

uu dt ^ uUf 

T ‘2uxUxxt 22ri Uxt T '^xxtt ^uutt — 0- 


(3.29) 


In this case L is nondegenerate, so that on © we can define the system of canonical coordinates 

Jg, =«(*-!) i=l, 2,3 

= *=1>2,3 

which reads 


qi 

= U 


<?2 

= Ux 


<?3 

— Uxx 


Pi 

— 22l/, Ux ^UxUxx Uxxt H“ SuUf 

P2 

= -6u^ 

“1“ Uxt “1“ ^UUxx 

P3 

= 6uux 

- Uf 


We will now solve the problem from the Hamiltonian point of view: starting from L and following Theorem 
1.1 one construct the canonical coordinates {pi,( 7 i} on ©: 


qi =u 

<?2 = Ux 
Q3 ~ '^xx 

Si 

Pi = 7- = TOu'^Ux - lAUxUxx - ^^UUxxx + 

Ollx 

P2 = 7- = ^‘^UUxx - U^^'> 

SUxx 


P3 = 


SI 

SxixX' 


— Uxxx 


and the Hamiltonian function 


Q = A- '^pi{q,)t. 
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By direct calculation one obtains 


35 

Q = 35u^ — lAOu^U^ - —u^ — + 70uu‘^Uxx — + 84u^Ua:'Ua;a:a: + 

- ISUxUxxUxxx - ^Ouulxx + + 6uUxxU^^'> - ^ + UxxxU^^^ 

and in canonical coordinates 

35 

Q = 359 ^ + 2SQqlql - —q^ - SSgr^ga + lOqiqlq^ - 21qlql + 

-ql- 6 qiq 2 Pi - &qlp 2 + ^qiqzP 2 - ^pI - ^^qlq 2 P 3 - 4 g 2 g 3 P 3 + PiPs + 4 gip|. 

The corresponding Hamiltonian system reads 

qi = 6 qiq 2 - Ps 

92 = 69 I - Sqiqs + P2 

93 = 70^192 + 4^293 - Pi - 891 P 3 

Pi = 21 Qql + 84091 gf - I4O9193 + 7O9I93 - 42 qiql - 692P1 + 893P2 - I4O9192P3 + 4p3 
P2 = 5609 i 92 - 7092 + I4O919293 - 691P1 - I292P2 - 709iP3 - 493^3 
P3 = -359i + 7 O 919 I - 429i9^ - 39 ^ + 891 P 2 - 492 P 3 

Rewriting this system in coordinates {pi, 9 i} one obtains exactly (3.29). 


3.A Appendix 


Proof of Lemma 3.1: we prove the Lemma in two parts: 


• firstly we prove the relation 


dk 

34> 


n 

= E 

j=i+m+l 


SI \ 


f = 0,..., n — m — 1. 


(a.l) 


For convenience we can explicitly rewrite eq. (3.8a), using (3.9): 


E / UJ^ \ UU"' 

) dt 



j^m 


E 


f dL \du^i'^ 


V J 9 ub) 


(d , 
ul ' + 


n—m r n 

■E E 

i—Q ^j—m+i 


/ dL 



1 dui^U 


(a.2) 
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Notice that the arguments of in the square brackets depend on u and its a;-derivatives upon the order m — 1 
and on ut and its x-derivatives upon the order n — m, so that the dependence on forn —m + l<i<m —1 

and on ul{\ for every j is only the explicit one. Analogously 


dx 


dA 


dA 


dx 

dA 

du 


dt 


E 


dA 


aub-i) 




Ut 


dA 


du. 


rn — l) 


dut 


-Utt 


E 

i— 

— m 

E 


i=l 

n—m 


dA 


dA 


.du.^ 

dA 


(b 


du(d J 


,b) 






(i) 


(a.3) 


The i-th step of (a.l) is obtained from the obvious identity 


g fdA\ 

du[f\dt) du\f\dx) 
Indeed, from (a.2) and (a.3), it follows that 


i = 1,..., n — m. 


d 



du^^^ 


dA 

duif Ua; 7 dx V duif / SuF ’ 

In particular, at the first step, i = n — m one obtains the basic relation 


dk 


du^'^'> 


au(”) 


\du^^^ J 





Substituting (a.6) into the further step of the recurrence, one finds 


dA 

f ^ ) 

I du^" , 

1 

f ^ ) 


d^in-m-2) - 

\Su(" ^Hx) J 


\Su("^ (x) J 

2 ) 


and so on. This gives relation (a.2). 


(a.4) 


(a.5) 


(a.6) 


• The second step is the proof of the relation 


dA 

EF 



which is a part of (a.l), indeed, for i > n — m, the partial derivative 
previous case, equation (a.7) follows from the identity 


i = n — m,... ,m — 1, iu-7) 


vanishes. Very much as in the 


d fdL\ d fdA\ 
dutt \dt ) dutt \ dx ) 
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Using (a.2) one can rewrite the left hand side as 


d 

dutt 



dut 


On the other hand, from (a.3), one has 


where 



dk 

^ dut ’ 

IT . \ 

(x) dut J 


We develop the right hand side, recalling (3.6b): 

dx \ dut ) dut dut 


obtaining 


d f dk 


dx V dutt 


= E 


d r TT \ r 'TT \ dT^^ f 'TT \ au(™) au(™) 

— 1 iw-T^ 1^=^ + dut \6u^i){x)) du^'^-^) dut 


(ix \(5'ut-^)(x) J du 

7T M du^ 


dx I (5M(’"+i)(a;) j 1 dut ^ 

^ ^ j—m+1 


df_ 


SI 


au(^) 


Y SI __ 

\(5m(-?)(x) ) ^ dx )\ dut 


- E 

j=rn+l 

±(- 


\( ' 

\ du^rn) 

auY)- 

\5u(d(^x)^ 

) auY-i) 

dut 


51 




/ dL \ dT^^ / (5/ \ dT"^^ du^’' 

1 7a\ I 1 T 77w 7 I 7^ TTTTTiT o 


E / \ ! 01 \ ou^" 

_ Vau^-’V dut \du^^T^)) dut 

3 —?7T. + 1 


Inserting in (a.9) and equating it to (a.8) we obtain 

Y au \ d f '51 \iauY) _ dk aY™) / YY \ aY™) aY”") 

\aMY)y dx \(5t6Y+i)(a;) y awt duT^~^^ \5u^^T^)) 


But the term 


is nonzero by definition, so that 


dk 



/ 51 \ 

/ a^Y-i) ’ 


which is the first step of the recurrence (a.7), and so on. 


• Finally, since (a.l) for i > n — m coincides with (a.7), it remains to prove it for i < n — m. 
relations can be obtained from (a.2) and (a.7) together with the identity 


a f^ 

a^W \dt)~ a„W V da; y ’ 


i = 1,..., n — m. 


(a.8) 


(a.9) 


These 
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Indeed, starting from the index i = n — m and using (3.9), one may write 


d 


du\ 


<^\ _ / 


dt 


dL 




dL 


d f dL 9m(") 


dL \ ^ f dL \ du^^') 

dt Vi 9 u(”) J ' ^ \du(d ) ' 


+ 


(a.lO) 


On the other hand 


d fdk 


dJ'^-^^\dxJ dx\dJ 


= L(. 


dA 


dA 


du. 

Performing the same steps as in the previous case, one obtains 


(n—m —1) 


+ 


dA 


du 


'm —1) * 


dA 


-if- 


dA 


SI 




+ E 


i JT \ 


which is the first recursive step of (a.l). 


Q.E.D. 


Proof of Lemma 3.2: We prove the Lemma by mean of the equivalence 


SI 

Su{x) 


dA 


dA 


duf^) 


which follows from the identity 

d fdL 
dut V dt 

Indeed, expanding, one has 

E- 

dut \dt ) dt dut du 

d ^ / dL i du^^f 
dt ^ / dut 


d /dA\ 


dut \dx J 


E 

du 


^ / du 


_d 

dt 


d i dA\ dA 
dx\dutt) dui"'~^^ dut J 


du) ^ ^ \dufd J du ’ 


where the last equality follows from the equivalence of (a.8) and (a.9). Expanding the right hand side 


d idL 
dut V dt 


On the other hand 


_d 

dx 


dA 
dt V duft 


di dA 

dt\dui"'~^^ )\ dut 


+ 


dA d fdu^'^^ 


du\ 


dt \ 


dut 


dL \ \ 

du ' ^ ^ 


j=m 


i dL \ du'^^'^ 
\duf^f) du 


(a.ll) 


d /dAS _ d / dA\ dA du^'^'^ dA du^^^ 
dut\dx) dx\dut)~^duf"^~d dut 9u 

d / dA \ dA du'^^f 
dx Vi9ut / duf^~d dut 


du 


dA d 
(’"-i) dt V 


dut 


W SI \ ^ \ 1 duf"^f dA d /dui 

\Suf'^''{x) J ^ \Suft)(x) J duf"^~^f du \ 






dA d fdu^^f 


(a.l2) 
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Comparing (a. 11) and (a. 12) one obtains 


jL 

dx 


dA 

dut 


-if- 


dA 


dt V dut 


dA 


dA 




du 


^ f dL \ du^^^ 


du 


51 


But Lemma 3.1 states that 


dA 

dut 


_±( dA 


dt \dutt 


§u(™){a 


51 

5ux{x) 


E 


51 \ 




du^^'i 

du 


- E 


so that 


E 

dx 


dA 

dut 


_±( dA 


dt V dutt 


= ±(- 


51 


r TT \ 

\(5u(l)(x) / du 


dx \ Sux{x) J dx I 

( TT 


E 


\Su^^'^{x) J \ du 


du^^-d 


+ E 

j=m+l 

=l(. 


d (du^^ 


y5u^d(^x) J dxi^ du 


51 


dx \5ux{x) 


E 

dx 


E 


f SI 


\Su^^^x) J\ du 




. E I (5uO(a;) / I i9u 

j=m+l ^ \ / / \ 


du 


Substituting in (a. 13) we get 

dA \ _±f dA 


^y(m-l) 


dt\dut^-^^ 


du^rn) 

dut 



du^d 


^ f \ du^"^ 

.E I (5uC)(a;)/I du du^'^~d 

j=m-\-l ^ \ / / \ 


du 


(a.l3) 


+ 


51 


[ \ (x) 


- E 

j—m-\-l 


du^^'l 


f 'll \ 

\i5u(l)(x)y J du 


All the terms cancels but 


dA 


dA 


du^^') 

du^rn-l) J dut 


-i( 


51 


dL\ 


du J dx \5ux{x) 


Q.E.D. 


Proof of Lemma 3.3: Relation (3.13) follows from the identities 

d (dL\ d (dA\ 

\lU ) ^ 


(a.l4) 
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and 


d f ^ 

du^f’ \ / du^t^ \dx ) 


Starting from (a. 14), we can write 


i = 1 ,..., TO — 1 . 


d fdL\ 
i9u(d \dt ) 


dt \9m(®) ) dt ) 9m(®) dt \i9u(d ) ’ 


j=m 


(a.l5) 


and 

d /dA\ _ d / dA \ dk dk dk d 

(9u(*) \dxy dx 9m(®) V du^'‘'' ) 

which give 


d 

• n 

r 1)1'' 

\ dud) 1 

^ dL Y 

dt 

E ( 

\dud) ^ 

1 dud) \ 

ydud) J 


dx \i9uh) y i9u(™-i) du^d \ du^d ) 


On the other hand, in (a. 15) 


d /dL\ dL 

Qu^d \dt) du^d 



du^^^ _d 
dud) dt 


dL 

du\^'^ 


and 

d /dk^^ _ d / dA \ dA dA du^"^) 

du[d \ dx) dx Vartl^V ^ dui'^~^'^ dud) 

which give 



dud) 

dud) 


'9A dA du^^) 

^^\du^d) ^^\du^d) du^^~d du['^~d dud) 


Performing the derivative w.r.t. t and substituting into (a. 16) one gets 


which gives 


d^ ( dL 


d d f dA 


dt^ \du[d 
d / dk \ 
dx \ dud) J 


dxdt \Qy^(d 

dk 


+ 


d dA 


dt du'^~d dt 


+ 


dk du)"^) 


dud-d du)"^~d dud) 


d \ dA du)^) 

dud) 

dk d f du^"^) 

dt \ dud) 


dL 
Kdu^t'^ 
dk 


jL 

dx 


\dud) J 


_£( dA 


dt \dui"'^J\ 


dA 


du)^-d dt 


du)^) 

dud) 


dA 


d dA 


dud-d dt Qyd~d \ 


The left hand side of (a. 17) is zero if t > n — to. 
If i < n — TO, 


dL 

n 

- V ( 

(m'' 

\ du)^) 

■ dk 

du£ 

2^ 1 

k—m-\-i 


^ du£ ~ 

ydut^^ 


+ 


dA 

dx \du^^ 


[a.n) 


(a.l8) 
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Indeed, using (a.l), and relation (3.6b), which we rewrite here below, 


dx V du["'’ ) du["^ awl*-!) 

one obtains 

aA V ( \ _dr ^ / 'IT 

a^-i) ^ ^ L)k+i au« . 



where the last identity follows from the fact that 

a^i^) 

if k — I < m. In (3.6), this implies 

g^{k+7n) g^{k+m+j) 

Finally, substituting (a. 18) into (a. 17) gives (3.13) 

' aA d / aA \ d^ / aA y _ 

_au(*-i) diV^yb-i); dt2 ivay(j-i) j 

_ ' aA d / aA \i d r aA d / aA \ d^ / aA y 

a^i’^-i) dt ) au^a dx aw^a dt \i9u|ay dt^ J 

Q.E.D. 


4. Applications to Painleve equations 


In this Section we study some applications of Theorem 2.1 and Theorem 3.1. to show how the finite 
dimensional Hamiltonian structure of Painleve equations comes from an infinite dimensional structure via 
the above procedure. 

4.1 PI as scaling reduction of KdV 

At the beginning we study the problem following the Hamiltonian scheme, then we will apply the 
framework of Theorem 3.1. 

We consider the KdV equation 

Ut = 6uUx - Uxxx- (4.1) 
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i.e. the t = tj-flow in the KdV hierarchy (1.2); it admits the nondegenerate scaling symmetry 


= J {u^ + ^ + + 6tu^)dx, 


(4.2) 


which depends on x,u,Ux,t. We note that L = [L^ + + 12tLp],where = f L_^dx = f and 

Ig = J L^dx = f “ dx are the first Hamiltonians of the KdV hierarchy. 

Theorem 2.1 states that the t-flow is Hamiltonian on the manifold © of the stationary points of the symmetry, 
i.e. 6 is the 2-dimensional manifold of the solutions of the Euler-Lagrange equation 


SI 

6u{x] 


= Ur.T. — in — 2x — 12tu = 0 . 


(4.3) 


It is invariant under the t-th flow and it naturally carries the system of canonical coordinates: 

q = u 

P ~ ^ ~ 


Notice that the identities 

^0, 

hold, where H is the generalized Legendre transform of L: 

SI 


(4.3) 


H = -L 


Su. 


-Ux- 


The first of identities (4.4) allows us to express the higher derivatives for m > 2 in terms of xd, p, q 
and with I = 1,..., to — 2 + 1. 

On & px + = 0, and the system (4.4) reduces to the canonical Hamiltonian system 

Px = 3(7^ + 2x + 12tq 
Qx =P, 


for the Hamiltonian function 


H = —L + ^- q^ — 2qx — 6tq^ 


(4.5) 


giving the reduced x-flow. This system is equivalent to the second order ODE in the variable q : 

q = —iq^ — 2x — 12tq. (4.6) 


The space © is the set of the stationary points of the scaling symmetry (4.2); this means that © carries a 
“natural” system of canonical coordinates {wi^TZi}, given by the self-similar function of u, i.e; combinations 
of u,x,t in the variable z{x,t) invariant w.r.t. the scaling. We will call them scaling coordinates. In this 


case 

f w = ^ + t 
I TT = 2p 


with z = X — 

In terms of the scaling coordinates the system reads 

{ dTT dS^ 

dz dw 

dw _ 
dz d-K 
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for the Hamiltonian 


f) = ^ - 8w^ - 4wz + 8t^ + Uz. 

8 

The system is equivalent to the ODE: 

w = + z, (4.7) 

that is exactly Painleve I. The Hamiltonian differs from the usual PI Hamiltonian for the terms in z, t 
that do not enter in the Hamiltonian system. 

We now construct the time dependent Hamiltonian function (—<5), that is the reduction on © of 

-0 = -A + p§, 
at 

where p, q are expressed in terms of m, u^, and A(a:, t, u, u^, Uxxx), calculated from 

dL dA 
dt dx 


has the form 


9 1 

A = 6^(4^^ - 2uUxx + ul) + 2x(3m^ - Uxx) + + “^^x - 8v?Uxx + - UxUxxx- 

By direct calculation one obtains 


9 

Q = 12^(2^^ H—^ — uuxx) H—^ ~ ‘^'^x + 2x{3u^ — v-xx)-, (4-^) 


This reduces on 6 to 


Q = — 6tq^ — 2xq) -\-2p — 2x^ 


(4.9) 


(4.10) 


Theorem 2.1 states that (—Q) is the Hamiltonian for the reduced t-flow,i.e., in terms of p and q 

(q=-2{6tp+l) = -§ 

\p= -m i3q^ + 2cc + 12tq) = 

Notice that system (4.10), written in terms of the scaling coordinates w and z, gives the same Painleve I. 

Remark: In this case the evolution equation is Hamiltonian and it can be written in the form 


where /^ = / L^dx with density 


i t \ T \ ^ 

«, = Wx),7.) = -^, 


T - 1,3 -L 

L, - u + 2 • 


On the other hand the scaling symmetry defines the stationary flow 

du 


ds 


— 12t'Ua; 4“ Ut 2 — 0. 


which is Hamiltonian: 


du r / N n ^ n 

— = {u{x),I} = —— = 0. 
as dx du 
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The s-flow and the t-fiow commute, but the Hamiltonian generating the scaling depends explicitly on 
the time t, so that the relation 

dL,. dl, 




= {Iuvh] + 


(«) 


= 0 


dt ' dt 

holds. Hence we have an alternative way to define the reduced Hamiltonian Q, following [BN]: 


d ^ dL dl d dll 
dx dt du dx du 


In this case the relation (4.2) follows as a consequence. 


System (4.10) i.e. the reduction of the f-flow on ©, can be obtained from the Lagrangian point of view; 
indeed one can consider the evolution equation (4.1) as the definition of Uxxx in terms of {u,Ux,Uxx,Ut), 
explicitly: 

^xxx — tSutlx Uf. 


Differentiating this relation in x one obtains 


= duuxx + Gul - Uxt 
m(5) = ISUxUxx + 36u^ - 6uUt - Uxxt 


which is a map from the ’’coordinates” 


u, Ux, Uxx, , 


into 


"U, Uxi 'Uxxj 'Ut: '^xt: 'f^xxtj 'Uttj ■ ■ ■ 


The Lagrangian L depends on u, Ux and hence its derivative ^ looks like 


dL _ dL dL dL 

dt dt du^* dux 


And there exist a functional A depending on u,Ux,Uxx,Ut such that 


dL dA dA dA dA 


dA 


JJL - J - Cl Cl a '^xx “f ^ ^ Cl ^xt' 

dt dx dx du dUx duxx dut 


dA 


Here, in terms of the new coordinates 

- 9 a 1 

A = 6t(4M® - 2uUxx + ul) + 2x(3u'^ - Uxx) + -j^U + -ulx + 2Ux - Su'^Uxx + UxUt, 


L = L, and 


dL „ 2 

a =®“ 

dL 

—— = 3u^ + 2x + 12tu 
du 

dL 

■ 5 — = 

UUx 
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= 


dA 
dx 

dA o o 

= 72tu"‘ — 12tUxx + 72xu + 18u'^ — Quux 
au 


dA 

dux 

dA 

dUxx 

dA 

dut 


— \2t'iix ~\~ lit “t“ 2 
= —V2tu + Uxx — 2x — 


The condition |^ = 0, that defines the submanifold 6, is equivalent to the condition 


dA 

dUxx 


= 0 


Hence we have an alternative definition of the space ©, and an alternative way to defines the canonical 
coordinates: 

q = u 


Theorem 3.1 states that the reduced t-fiow is Lagrangian ,with Lagrangian A, in this case it is easy to verify 
it, indeed, on 6, 


= —12tu + Uxx — 2x — 3m^ = 0 

= 12tUx + Ut + 2 = 0 

M - 72tUxx + 12xu + 18u^ - Quuxx - Uxt = 0, 

where the first equation is the definition of the submanifold © itself, the other two reproduces (4.9), indeed 
they can be rewritten as 

( Ut = —12tUx — 2 
1 Uxt = —12t{8u‘^ + 2x + 12tu) 


4.2 PII as scaling reduction of mKdV 

One can repeat the same procedure as in section 4.1 starting from the mKdV equation 

Ut = Qu^Ux — Uxxx- (4-11) 


It admits the nondegenerate scaling symmetry 

f , A 9 ^ u^x. , 

I = J {-t{u +Ux)-\ ^)da;, 

2 

which depends on x, u, Ux, t. We notice that L = 3tL^ + 

Here © is the 2-dimensional manifold of the solutions of the Euler-Lagrange equation 

Si 1 

. / N = Uxx - + ux) = 0. 

du[x) 3t 


(4.12) 


(4.13) 
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It naturally carries the system of canonical coordinates: 


q = u 

P=m 


— = 3tUr- 


As in the previous case we read the Euler-Lagrange equation as a reduced x-flow with Hamiltonian 




(4.14) 


where H is the generalized Legendre transform of L: 


SI 

H = —L + -—Ux = —L + 3tw 

OUr, 


2 

X • 


The system is equivalent to the second order ODE in the variable q : 

s 1 

qxx = 2g + —(jx. 

The scaling coordinates are now 



in the variable 2 ; = —and the system transforms into 
(3i)3 




( dTT dS} 

) dz dw 

1 d^ _ dS^ 




\ dz diT 


for the Hamiltonian 


! 2 4\ 1 2 


The system is equivalent to the ODE: 

w 

= + ZW. 

(4.15) 


that is exactly Painleve 11. 


We now construct the time dependent Hamiltonian function (—Q), that is the reduction on © of 

dq 


—Q — ~A + p 


dt 


where 


A = - Qu^Uxx + - 3UxUxxx) + + 2^1 ~ UUxx) + UUx- 


By direct calculation one obtains 


which on © reduces to 


13 1 

Q = - u^Uxx + + ^( 2 “^ “ 


PS ^(P 4 n ,1 1 2 2 


and is the Hamiltonian for the reduced t-flow . In fact 


p, — —JL — 

y 3t qF dp 

P, - p. - SeI - On^T - ^ 

.P~3t 3i , 


(4.16) 

(4.17) 


(4.18) 
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Notice that also the system (4.18), written in w and z, gives Painleve II. 
Remark: The evolution equation is Hamiltonian and can be written in the form 


where /^ = / with density 


I ( ^ T \ ^ 


i. = -(«* + »!) 


On the other hand the scaling symmetry defines the Hamiltonian stationary flow 


du d SI. 

— = XUx + AtUt + U= - -p 

as ax ou 


(s) 


= 0 . 


We now deduce system (4.18) from the Lagrangian point of view, reading the evolution equation (4.11) as 
the definition of Uxxx in terms of {u,Ux,Uxx,Ut), explicitly: 

Uxxx = Gu^Ux - Ut- 


Differentiating this relation in x one obtains 

J = Gu^Uxx + 12ut6^ — Uxt 

\ = iGuUxUxx + iGu^Ux + 12u^ — Gv?Ut — Uxxt 


which is a map from the ’’coordinates” 


u, Ux, Uxx, , 


into 


W, Uxj Uxxj Uf, Uxt: Uxxt: ^tt: • • ■ 


Here 

A = t{Gu^ - Gu^Uxx + -julx 

and 


3 / 4 , 2 \ 

-{u +Ux)-UUxx 

36tu^ — IStu^Uxx + Gxu^ — XUxx + Ux 
XUx + M + 3tut 
3tuxx — Gtu^ — ux 
3tux 


dA 

dx 

dA 

du 

dA 

dux 

dA 

dUxx 

dA 

dut 


3uxUxxx T 3uxUtJ 4” ^( 2 ^ 4“ 2 ^^ uuxx') 4“ 


UUx 


dL 8 4 2 n 

» = 2 <" +“-) 


du 

dL 

dux 


= Gtu^ + ux 
= 3tux 
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The condition = 0, that defines the submanifold 6 , is equivalent to the condition 

dA 


Out 


= StUrx — QtU — UX = 0 


Hence we have an alternative definition of the space ©, and an alternative way to defines the canonical 
coordinates: 

q = u 

P=§^= 


On 6 

f ^ ^tUxx - - UX = 0 

i ^ = xux + 3tut + u = 0 

~ 4i (^) = -IStu'^Uxx + + Gxu^ - xuxx - 2 mx - StUxt = 0 , 

where the first defines the submanifold, the second one gives the motion of u and the third the motion of 
Ux, hence one can rewrite them as 

f itUt = —XUx — u 
I “itUxt =Ux- 2u^x - 

which coincides with (4.18). 


4.3 PHI as scaling reduction of Sine-Gordon 

A particular case of Painleve III equation can be obtained as reduction of the Sine-Gordon equation 

(4.19) 


^t = v=^ 
vt = Uxx - sm u = 


via the scaling 


^ = ^y + tUx =^-^=0 

dv 


SI, 


SI 

= x[Uxx - sm U) + tVx +Ux = — 5 ) 7 - = 0, 


where j L-^^dx, I,, = f Ldx, with the Hamiltonians 


(4.20) 


and 


w.r.t. the Poisson bracket 


+Ux:)-cos 


L = + Ux) — XCOS U + tvUx = xL^ + tvUx 


{F.G} = Ji 


Sf 6g 6f 6g 


du{x) 5v{x) Sv(x) Su(x) 


) dx. 


The scaling reduction equation means 


which defines the submanifold ©: 


with the canonical coordinates 


K) ^ K) 0 

Su(x) Sv(x) 


t ’ 


+ 2 _ 2 

p = XUx tv = ^ V 

q = u 


(4.21) 
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in 6. 

The equation defining © can be written as an Hamiltonian system in canonical form describing the 
reduced x-fiow: 


= -f 

(9)x = ^ 


where 


In terms of the scaling coordinates 


TJ ^ P , 

U = — — -^ + X cos q. 


2 x"^ — 

{ 


(4.22) 


w = q 

-2z 

7T = -rP 


2 _.2 

in the variable 2 = ^ ^ , the Hamiltonian system transform into 


d7T dSj 

dz dw 

dw dS) 

dz Ott 


for the Hamiltonian 




The system is equivalent to Painleve HI: 


1 TT^ 

'iV 


— cos 


w. 


2zw + 2w — sin w = 0. 


Let us now construct the time dependent Hamiltonian function {—Q), that is the reduction on © of 

-<3 = -A + P§. 

at 


where 


By direct calculation one obtains 


which on © reduces to 


A = xUxV + + v?x) ~ cosu). 




1 

Q = -t inT2 -12 - COS q). 


'2 

This is the Hamiltonian for the reduced t-flow . In fact 


q = V = --Ux = —t ,2^ 2 = — a 

^ X t^ — x^ op 

p = V + XVX + tvt = t sin q = 


(4.23) 


Note that also the system (4.23), written in w and z, gives Painleve HI. 

Remark : We now deduce system (4.23) from the Lagrangian point of view, reading the evolution equation 
(4.19) as the definition of v in terms of Ut, explicitly: 


( V = Ut 

\ Uxx = ft + sin u = utt + sin u 
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Differentiating this relation in x one obtains 


Vx — y^xt 

Vxx =Vtt-V COS U = Uttt - Ut COS U 


^XXX — Uttt + U^cosu 


which is a map from the ’’coordinates” 


into 


"U, n, Ux-, Vxi Uxx: Uxxj ■ 


"U, Uxj Utj Uxti Uttj Uxttj • 


Here 


and 


which give 


and 


A = xUxUt + — + u^) — cos u) 


L = tUxUt + x(-{Uf + Ux) — cosu) 


w 

du 

dL 

dux 

dL 

dut 

dA 

dx 

dk 

du 

dk 

dux 

dk 

dut 


2 - cosu 


= — sm u 


= tUt + XUx 


= tUx + XUt 


2 - cosu 


= —tsinu 


= XUt + tUx 


= XUx + tUt 


The condition 4^ = 0, that defines the submanifold 6, is equivalent to the condition 


dA 

duo, 


= XUt + tUx = 0 


Hence we have an alternative definition of the space ©, and an alternative way to defines the canonical 
coordinates: 

q = u 
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On 6 


I ^=xut + tu^ = 0 

I - HM) = sinu - xu^t -ut- tutt = 0, 
where the first defines the submanifold, and the second 


—t sinu — xUxt — Ut — tutt = 0 


coincides with (4.23). 


5. Self—similar solutions of n—waves equation and Hamiltonian MPDEs 


5.1 n—w8cves equations and their symmetries 


Let us consider the equation 

Ut - Vx - [u,v] = 0, (5.1) 

where 

u=['),a] v=['y,b] a = diag {a}, b = diag {b^, (5-2) 

and 7 is a function of x, t. 

Following [DS] it is possible to rewrite (5.1) as an infinite dimensional Hamiltonian system on the 
space 9Jt of functions of x with values in Mat{n, C) with the inner product 


(■u,u) 


Tr {u{x)v{x))dx. 


On the space of functionals 


F = 


ft 

}[X,U,Ux, . — u 


) dx 


one can define G 911 by 

^F(u + ew)L^o = (V„F’,w) 
and the Poisson structure P with the Poisson bracket 


{F,G}iu) 




(5.3) 


The n-waves equation (5.1) is a Hamiltonian system w.r.t. this Poisson structure: 


Ut = Pdli 




A 

dx 


■ u\ = \—v 


’ dx 


(5.4) 


where 


that in components of 7 gives 


h 


Lidx 


1 

2 


Tr (uv)dx 


h 


= [ — T — u^kUki)dx= ( 

J j J ■ j 


(5.5) 
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T 

For {n = 3, u = —u) one can reduce to a particular case of P VI equation (see [D], where the 
Hamiltonian structure for this particular case of P VI is derived from the Hamiltonian structure of the 
n-waves equation.), imposing the scaling 


It admits the Hamiltonian form 


where 



du 

— =tUt+ XUx + u = 0. 
ds 


du ^ d 1 „ 


Tr {tuv + xu^)dx = / {t -— 

J ^ “ ak- ai 

I k 


2 )uikUi^idx 


and 

or, in terms of 7 : 


V„/(s) = —tv — XU, 


hs) = /EEk^* - h){a^ - ak)t + (a* - afc)^x] 7 ,,i^ 7 ^, 


dx 


(5.6) 


(5.7) 


(5.8) 

(5.9) 


We emphasize the fact that the t-flow and the s-flow commute, so that 




dl, 


(o) 


dt 


= 0 . 


By substituting: 

/ d 1 

[Tr(V„/i[V„/(s), — + m]) + -Tr{uv)]dx = 0. 

Then there exists a function (x, t, u, v) such that 


d uv 

Tr{—v[—tv — XU, ——h mJ + —) 



(t) ■ 


By direct calculation (see Appendix 5.A) we obtain 


Q(,t) 


-Tr{xuv + tv^) 


I E[(a^ - -b,)x + {b. - 6j"t]7y7,. 

i,3 


As in the previous examples, is the Hamiltonian for the reduced f-flow. We now describe this 

flow. 

We start by rewriting the system 


in terms of 7 , i.e. we solve 


under the condition 


This gives 


{ut-Vx- [ m , w ] = 0 
\ tUt + XUx + u = 0 


[it, a] = [ 7 „, 5 ] + [[ 7 , 0 ], [ 7 , 6 ]] 
t I 


[ 7 t, ax + tb] + [ 7 ,5] = [[ 7 , ax], [ 7 , 6 ]] 


(5.10) 
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but, because of the commutativity of b with itself, 


dt 


[ 7 , ax + tb] = [[ 7 , ax + bt], [ 7 ,5]]. 


(5.11) 


Then we identify &s with the space of matrices 

q = [ 7 , ax + bt] = ux + vt, 

so that u = ad, ,_,ad“^q, v = ad, ,_,adr^q, or, in terms of the matrix elements: 

dij = [(a^ -cii)x+ (bj - 

On © the equation (5.11) has the Lax form 


with the Hamiltonian function 


This coincides with . 


qt = [q,v] = [g,ad(^^^^,^,ad/q] 


= ^Tr(gu) = ]^Tt{xuv + tv^). 


(5.12) 


One may change the role of x and t. This means that one considers the system (5.1) on the space of functions 
v{t): 

Vx = +u], 

where do = f Hodt = f Ti(uv) dt and one integrates in the variable t. The scaling (5.8) can be read as 
an Hamiltonian equation 


where 


and 


= 



o' 

II 

+ 

(5.13) 


Tr (xuv + tv'^) dt 

(5.14) 

V J(.) = ■ 

— tv — XU 

(5.15) 


Commutativity of the flows is equivalent to 


dl, 


(7 / I-Hi = 0- 

O0O(,)l 


in our case 


/ T 1 

[Tr(V„/o[V„J(,), j^+v]) + ^Triuv)]dt = 0. 
Then there exists a function {x, t, u, v) such that 


, r d . uv. d 

Tr(-«[-tu-xr.,-+u] + -) = --Q^^,. 


By direct calculation (see Appendix 5.A) we obtain 


1 . 


Q(:c) = ^Tr(tuv + XV?), 


(5.16) 
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in components: 


Q 


(x) 


I ^[(a, - aj( 6 , -b.)t + {a. - 

i,3 


Now we study the a;-flow on the reduced manifold defined by the scaling equation : 
the system (5.10) gives 

{ut-Vx- [ m , w ] =0 
\ tVt + XVx + u =0. 

In terms of 7 this becomes 

— [ 7 , aa; + t 6 ] = [[ 7 , ax + bt] ,[ 7 , 0 ]], 

that is a Lax equation on ©gi 

qx = [q,u] = [ 9 ,ad(^^^^,^,ad;^q] 

with Hamiltonian function 

= ^Tr(gu) = iTr(xu 2 + tuv). 

This coincides with ■ 

In fact one can rewrite the scaling as a zero-curvature equation in two ways: 


du 

ds 


qx + [u, g] = 0 


and 


du 

ds 


qt + [u, q] = 0 . 


Therefore one may rewrite them in terms of q as 


9 - = ld,ad^^l^^^^ad^q] 


and 


qt = [q,ad^2.^^^^ad,q]. 


5.2 Commuting time—dependent Hamiltonian flows on so(n) 
We can do exactly the same using the coordinates 

ti = xui + tbi, 

and the corresponding derivatives with 


d \—A d 

= z.« 


dx 


dt. 


and 


The starting equation is now 


d d 


dt 


dt. 


- 4-Uk - [Mj, Mfc] = 0 

dtf- dt. 


(5.17) 


(5.18) 

(5.19) 


(5.20) 

(5.21) 


(5.22) 
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where 


M, = [ 7 , E,] {Ui)kl = 
and {Ei)ki = SikSki- We impose the scaling 


d ^—A d 

-f'^k = 2_^U-^Uk + Uk = 0 


(5.23) 


For every k one can define, on the space of functionals 


^ /■ ., du d^u. , 


—F{uk + ew)\^^„ = (V„^F,w), 


a Poisson structure P with the Poisson bracket 


{F, G}iuk) = (V„^ F, [V„^ G, — + uk]) 


The n-waves equation (5.22) is Hamiltonian w.r.t. the Poisson structure P in 5^^,: 

= [-u,, — +Uk], 


where 


/i = y Lidt = -i y Tr (uiUk) dtk- 


(5.24) 


(5.25) 


On we can reduce to P VI equation imposing the scaling (5.23), which admits the Hamiltonian 


, Ufc — [V„^/(s), + Uk] — [ + Uk] — 0 


(5.26) 


where 


^{s) — f Edtf^ — f Tr tjUjUkdtfi^ — 


tjUjUk + tkul)dtk 


(5.27) 


The commutativity of the flows is equivalent to 




m our case 


y [Tr(V„^/i[V„^/(s), -^+Uk]) + ^Tr{uiUk)]dtk = - j dkQ^.^dtk- 


By direct calculation (following the scheme in Appendix 5.A) we obtain 


‘^(i) = FTr^tjUjM* = 


(5.28) 


The scaling equation defines the submanifold ©«. One can consider on ©^ the system of coordinates 
given by the matrix elements of q: 

q = b,'^tjEj] = h,U], 
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where U is the diagonal matrix diag(ti,..., tn); explicitly 


= {tj (5.29) 

As in the previous cases, Q(i) is Hamiltonian for the ti-flow on the reduced manifold. 

Starting now from i k, we can reduce on the same submanifold © and construct the Hamiltonian 
function . 

Indeed, the scaling (5.23) for every k produces on ©^ the Lax equation 

qk = [q,uk], (5.30) 

with Hamiltonian functions 

Hk = ^Tr(gufc) = ^TiY^tjUjUk = ^ 


These coincide with the 


T 

In the case q = —q eqs. 
differential operator 


constructed above. Observing that 7 = ad~^q one can rewrite 
Uk = ad^j^adO^q. 

(5.30) are the Monodromy Preserving Deformation equations for the linear 




that give Painleve VI, for n = 3, and the higher-order analogues, for n > 3. 


Remark: The first integrals of the MPDE (5.30) are given by the monodromy data of the operator A. The 
Poisson bracket on the space of the monodromy data has been computed in [Ug]. 


5.A Appendix 

Here we present the explicit calculations giving rise to equation (5.16): let us consider the following 
explicit expressions: 


A = y Tr {uv)dx 

Vlt = -V 

/(s) = — - J Tr {xu^ +tuv)da 


V/(s) = —tv — XU 


{It,I(s)} = (-V,tUt +XUx +u) = 

= — y Tr {tvvx + tv[u, w] + xvux + uv)da 

dl, s If 

= It = --jTi (uv). 

1 


dt 

dl 


In (a.l) the relations 


dl f 1 

{/(, /(s)}-^ = - / Tr {tvVx + tv[u, v] + xvUx + -uv)dx 


Tr V [u, v] = 0 

Tr (vvx) = i-^Tr (u^) 

^ ’ 2dx ^ ' 

Tr (xvUx H— uv) = -—Tr (xuv) 
2 2 dx 


(a.l) 
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hold. In fact, in terms of 7 ^ one can write 


which implies 

Then: 


Tr (xvux) = EE x{bk - h){a^ - ak)lik{l3,)ki = 

i k 

= EE x{h - bk){ak - ai){'yx)kilik = 

i k 

= Tr (xvxu), 

^-Tr (xuv) = 2Tr (xvUx) + Tr (uv). 
ax 

= ^^r{xuv + tv^)dx. 
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